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ABSTRACT 

Context. Mode identification is a crucial step to comparing observed frequencies with theoretical ones. However, it has proven to 
be particularly difficult in rapidly rotating stars. An important reason for this is the lack of simple frequency patterns such as those 
present in solar-type pulsators. This problem is further aggravated in 6 Scuti stars by their particularly rich frequency spectra. 
Aims. As a first step to obtaining further observational constraints towards mode identification in rapid rotators, we aim to accurately 
calculate mode visibilities and amplitude ratios while fully taking into account the effects of rotation. 

Methods. We derive the relevant equations for calculating mode visibilities in different photometric bands while fully taking into 
account the geometric distortion from both the centrifugal deformation and the pulsation modes, the variations in effective gravity, 
and an approximate treatment of the temperature variations, given the adiabatic nature of the pulsation modes. These equations are 
then applied to 2D oscillation modes, calculated using the TOP code (Two-dimension Oscillation Program), in fully distorted 2D 
models based on the SCF (Self-Consistent Field) method. The specific intensities come from a grid of Kurucz atmospheres, thereby 
taking into account limb and gravity darkening. 

Results. We obtain mode visibilities and amplitude ratios for 2 M models with rotation rates ranging from to 80% of the criti- 
cal rotation rate. Based on these calculations, we confirm a number of results from earlier studies, such as the increased visibility 
of numerous chaotic modes at sufficient rotation rates, the simpler frequency spectra with dominant island modes for pole-on con- 
figurations, or the dependence of amplitude ratios on inclination and azimuthal order in rotating stars. In addition, we explain how 
the geometric shape of the star leads to a smaller contrast between pole-on and equator-on visibilities of equatorially-focused island 
modes. We also show that modes with similar (I, \m\) values frequently have similar amplitude ratios, even in the most rapidly rotating 
models. 

Key words, stars: oscillations (including pulsations) - stars: rotation 



1. Introduction 



The space missions CoRoT (Bag lin et al.| |2009] |Auvergne et al.| |2009) and Kepler ( |Borucki et al.| [2009) are revealing very rich 
pulsational spectra in rapidly rotating 6 Scuti stars. For instance, several hundred individual frequencies have been found in HD 
50844 a nd HD 181555, obse rved by CoRoT ( |Poretti et aT]|2009| Michel, private communication), and V2367 Cyg, observed by 
Kepler (Balo na~et al.| |2fJl2] l. It is becoming increasingly clear that interpreting these spectra will not be a straightforward task 
and that theory is lagging behind observations. A crucial first step in interpreting this data is correctly identifying the pulsation 
modes, i.e. finding the correct correspondence between theoretically calculated modes and observed pulsations. Recently, [Reese 
et al. (2009b) proposed a way to identify acoustic pulsation modes in rapidly rotating stars based on an asymptotic formula which 
describes the frequencies of low degree modes (see Pasek et al. 2012 and references therein). Nonetheless this method runs into 
trouble if chaotic modes are present in the pulsation spectra, which is expected based on the visibility calculations in Lignieres & 
|Geor geot (2009}. Furthermore, the pulsation modes in 6 Scuti stars tend to be of low radial order, and may therefore be too far from 
the asymptotic regime. [Lignieres et al. ( 2010| l have worked on using the cross-correlation of pulsation spectra. Although it doesn't 
yield individual mode identifications, it may provide a way of obtaining the rotation rate and/or the large frequency separation, and 



explaining recent observations of recurring frequency spacings in rapid rotators (Garcia Hernandez et al. 2009 Mantegazza et al. 
2012). Nonetheless, the need remains for methods capable of identifying individual pulsation modes 



Two particularly promising methods for identifying pulsation modes are multi-colour photometric and spectroscopic mode 
identification. The first approach consists in measuring the amplitudes and phases of a given pulsation mode in different photometric 
bands, calculating the ratios of the different amplitudes and/or the phase differences, and comparing these to theoretical predictions. 
The second approach exploits the Doppler shifts caused by the velocity field from the pulsation mode and how it affects observed 
absorption lines. These methods have been successfully applied to slowly rotating stars (e.g. |De Ridder et al.[ |2004[ |Zima et aL] 
2006 Briquet et al. 2QQ7\, but more work is needed before they are applied to rapid rotators. In the present paper, we will focus 



on mode visibilities in different photometric bands as a first step to multi-colour photometric mode identification, and postpone 
spectroscopic mode identification to a later paper. 
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Few studies have dealt with multi-colour photometric mode signatures in rapidly rotating stars, and those that do generally 
approximate the effects of rotation on the pulsation modes. For instance, Daszyhska-Daszki ewicz et aL| (2002 |2007| l and Townsend 



(2003a) used either the perturbative approach or the traditional approximation to calculate their pulsation modes, and in some 
cases included the effects of avoided crossings. Their calculations of mode visibilities included stellar surface distortion from the 
pulsation modes and the Lagrangian perturbations to both the effective temperature and gravity. An interesting result from these 
studies is that contrarily to non-rotating stars, amplitude ratios depend both on m, the azimuthal order of the pulsation mode, and 
i, the inclination of the star. More recently, Lignieres et al. (2006) and Lignieres & Georgeot (2009 ) calculated geometrical disk- 
integration factors of acoustic modes in deformed polytropic models by integrating the temperature fluctuations over the visible 
disk. The effects of rotation were fully taken into account in the pulsation modes, thanks to the 2D numerical approach, but non- 
adiabatic effects were neglected, thereby making the fluctuations of the effective temperature inaccessible. Furthermore, gravity and 
limb darkening, the Lagrangian perturbations to the effective gravity, and surface distortion caused by the modes were not taken 
into account. Nonetheless, important first results were obtained through these articles, namely, that chaotic modes are more visible 
than their non-rotating counterparts due to irregular latitudinal node placement and may thus be detected, island modes are the most 
visible modes in a pole-on configuration, leading to a regular frequency pattern (Lignieres & Georgeot 2009), and signatures of 



the large frequency separation and/or rotation rate can show up in the autocorrelation function of the frequency spectrum (Lignieres 

|etaL]|20T0l ) 

In order to obtain more realistic multi-colour photometric mode visibilities in rapidly rotating stars, we derive a new set of 
equations which take into account the Lagrangian variations to the effective temperature and gravity, as well as the surface distortions 
induced both by the centrifugal deformation of the equilibrium model and by the pulsation modes. These are applied to adiabatic 
acoustic modes calculated by the Two-dimension Oscillation Program (TOP, Rees e et aL] 2006 , 2009a i using rapidly rotating zero- 
age main-sequence (ZAMS) models based on the Self-Consistent Field (SCF) method ( Jackson et al.| |2005| |MacGregor et al.| 
2007 ). The emergent intensities are calculated from Kurucz atmospheres, taking into account the latitudinal dependence of the 
equilibrium effective temperature and gravity, thereby including gravity and limb darkening. The main weakness in the present 
study is the adiabatic approximation, which makes the Lagrangian fluctuations of the effective temperature inaccessible. As was 
previously done in Lignieres & Georgeot (2009 ), we approximate these by the Lagrangian temperature variations. The following 
section describes the pulsation calculations, with an emphasis on the improvements and differences with the calculations done in 



Reese et al. (2009a). This section is followed by a derivation of the relevant equations for calculating mode visibilities in rapidly 
4]then describes various effects of rotation on visibilities in a single band - the CoRoT photometric band. This 
which deals with amplitude ratios in the Geneva photometric system. The paper ends with a short conclusion. 



rotating stars. Section 
is followed by Sect.|5 



2. Pulsation calculations 

In what follows, we review the methods used for obtaining the models and associated pulsations that serve as inputs to the visibility 
calculations. These closely follow the approach used in |Reese et al.| ( p009a[ ) but include a number of improvements as described 
below. 



2. 1 . Equilibrium models 

The equilibrium models are calculated via the SCF method (Jackson et al. 2005 MacGregor et al. 2007| ). This method is an iterative 
procedure which alternates between solving Poisson's equation and the equations of mass, momentum and energy conservation 
before converging onto a 2D centrifugally deformed stellar model. These models are chemically homogeneous ZAMS models with 
a cylindrical rotation profile, although throughout the rest of the article, we will work with uniformly rotating SCF models, even if 
the formulas in the visibility calculations are established for general (non-cylindrical) rotation profiles. Given the rotation profile, 
the structure of the model is barotropic, i.e. all thermodynamic quantities remain constant on isopotentials, which are calculated 
from the sum of the gravitational and centrifugal potentials. Finally, we wish to make the distinction between the critical rotation 
rate, Oc, and the Keplerian break-up rotation rate, Ok: 




Ok 




(1) 



Although very similar, the former uses the true gravity (excluding the centrifugal force) at the equator, g eq , based on the actual 
distribution of matter, to calculate the break-up rotation rate, whereas the latter uses its Keplerian approximation, ^M, which amounts 



to assuming spherical symmetry for the distribution of matter. As pointed out in Roxburgh (2004), the Keplerian approximation 
slightly underestimates the true gravity, so that fic > Ok- Table [T] gives the relative differences between these two quantities for 
different rotation rates, calculated two different ways. The first method is based on global quantities provided with the models, 
whereas the second involves recalculating the gravitational potential from the density distribution and using this to calculate the 
equatorial gravity. A comparison of columns two and three, and five and six gives an idea of the uncertainty on these values. We 
also note that the theoretical value of (fic _ ^k) /^k at = 0.00 is 0, since the star is spherically symmetric and the Keplerian 
approximation is exact. 

Before being used in the pulsation calculations, the models need to be interpolated onto a new grid, and a number of supplemen- 
tary equilibrium quantities have to be derived, including a variety of geometric terms as well as gradients of different equilibrium 
quantities. Since Reese et al. ( 2009a| ), a number of improvements have been incorporated into these procedures. For instance, the 
stellar models are now interpolated onto a non-uniform radial grid which becomes dense near the stellar surface. This allows the 
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Table 1. Relative differences, (Qc - Qk)/^k, for selected rotation rates, calculated with two different methods (see text for details). 
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pulsation code to correctly resolve the rapid spatial variations of acoustic modes near the surface, resulting from the decrease in 
sound velocity. Instead of interpolating the density and pressure directly, their logarithm is interpolated. This leads to more accurate 
and consistent values near the surface, where these quantities are several orders of magnitude smaller than in the centre, and ensures 
they remain positive. The effective gravity is calculated via Poisson's equation. This avoids taking the ratio of the pressure gradient 
divided by the density, both of which are small quantities subject to relatively large uncertainties. Furthermore, the equipotentials, 
and hence the geometric structure of the star, are recalculated using the solution from Poisson's equation thus removing some nu- 



merical inaccuracies in the original models. The £ derivative of equilibrium quantities is now correctly calculated. In Reese et al. 
(2009a ), the derivative was mistakenly calculated with respect to r rather than Although this changed the quantitative results, the 



qualitative conclusions from Reese et aL"l(|2009a} remain unaltered. Finally, the Fi profile is not derived from Eq. (16) of Jackson 



et al. (2005 1, but rather from the equation of state, which is based on the formula of Eggleton et al.|(|1973[l. Furthermore, the quantity 



T2, which intervenes in the visibility calculations described below, is also calculated via the equation of state. 



2.2. Pulsation equations 

A new set of variables is used in the pulsation equations: 

€. % 'A % (2) 
Pq Po 

where £ is the Lagrangian displacement, Sp/Po the Lagrangian pressure perturbation divided by the equilibrium pressure profile, 
£p/po the Lagrangian density perturbation divided by the equilibrium density profile, and *F the Eulerian perturbation to the gravita- 
tional potential. Throughout this article, the subscript "0" denotes equilibrium quantities. We assume that the time and (f> dependence 
of these variables takes on the form e l (< Jt+m <!>\ where m is the azimuthal order. As such, we use what could be called the "retrograde 
convention", i.e. modes with positive azimuthal orders, m, are retrograde. 
Based on these variables, the continuity equation becomes 

dp 

= - + V ■ £ (3) 

Po 

and Poisson's equation is 

= A¥-4*gL£-£-VpoJ, (4) 

where G is the gravitational constant. Euler's equation takes some more manipulations (see, for example, Eq. A. 3 of |Reese et al.| 
[2009il l: 

-[oj + mQ] 2 £ + 2i [u + mfi] £1 x £ + &sd s fo 2 )e s = - + - V*F 

v ' Po Po 

= --v(po S 1 ~ € ■ VPo) + (po^ - f ■ V J ^ - V* 



Po \ Po j \ Po I Pq 

= P °y( 6p ) 6pVP ° i V ^- VP ») i s P* p o (£-Vpo)VA) v>1 , 
Po \Po) Po Po po Po Po pi 

. ft t m J » - f) ,, + v tf . - vr J 1 - -<* ■ *»^ F « } , ( 5 , 

Po \Po) \po Pol I pi 1 



where Q. is the rotation profile, s the distance to the rotation axis, and g e g = -V^o + sQre s = VPo/po the effective gravity. The term 
in curly brackets cancels out because VPq is parallel to Vpo in a barotropic stellar structure. Finally, the adiabatic relation takes on 
the following very simple form: 

* = 1* (6) 

po n p 

This last equation is then used to eliminate 6p/po in favour of 6p/Po throughout the differential system and thus to reduce the size 
of the problem compared to what is obtained in|Reese et al. (2009a). Explicit expressions for Eqs. ([3]), (|4]) and |5]), using spheroidal 



coordinates (see Sect. |2.4| ), are given in App. |A.l| 
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Besides reducing the computational cost, using the variables in Eq. Q allows us to obtain a much cleaner derivation of Sp/Po, 
and hence 5T/To, near the surface. Indeed, if one were to calculate this quantity from the Eulerian pressure perturbation, they would 
apply the following relation: 

dp = p + £- VPp 
Po Po 

Near the surface, this involves the sum of two nearly opposite terms, divided by a small quantity, thus leading to poor numerical 
results. 



2.3. Non-dimensionalisation 

Contrarily to |Reese et al . (2009a ), we non-dimensionalise the SCF models in a more classical way. The following reference quanti- 
ties are used as units of length, density and pressure: 



Pref 



M 

^eq 



P r 



■et 



GM 2 



(7) 



where R eq is the equatorial radius and M the mass. These lead to the following time scale: 



f ref 



GM 



1 



(8) 



where Qk is the Keplerian break-up rotation rate. Hence, the non-dimensional frequencies are directly oi/Q.^. With this non- 
dimensionalisation, the preceding pulsation equations remain unchanged except for Poisson's equation, which becomes 



= AT-47r|p(r^-£-Vpo 



(9) 



2.4. Spheroidal geometry 



As was done in |Lignieres et al. ( 2006 ) and Reese et al. ( 2006| l, a surface-fitting coordinate system, (£, 6, <p), based on Bonazzola 
et al. ( 1998}, is introduced. This system is related to the usual spherical coordinates, (r, 6, (p), via the relation 



r(£, 6) = (1 - e)t + 5 ^ - (R s (6) -l+s). 



(10) 



where R s (6) corresponds to the surface and s = 1 - R p /R eq is a measure of the oblateness, R p being the polar radius, and £ is 
comprised between and 1. As can be seen, £ = 1 corresponds to the stellar surface. A second domain, with a spherical outer 
boundary, is added around the star so as to simplify the boundary condition on the gravitational potential: 



■% 2 + \2^-A){R s {B)-\-e), 



(11) 



where £ € [1,2]. For £ = 1, Eq. (Hi coincides with the stellar surface, whereas for £ — 2, it yields a sphere of radius 2 (or 2R eq in 



dimensional form). For conciseness, we will use the subscripts "f" and "ff" to denote derivatives of r with respect to these variables. 
For example, r% = dr/d£ and r (g = d 2 r/d£d9. 



2.5. Boundary conditions 

The pulsation equations are supplemented by a number of boundary equations. Regularity of the solutions is imposed in the centre. 
The perturbation to the gravitational potential is made to go to zero at an infinite distance from the star. This condition is imposed 
by extending *F into the second domain and matching it to a vacuum potential on the second domain's outer spherical boundary as 



described in Reese et al. ( 2006| >. When extending *F into the second domain, both it and its gradient need to be kept continuous across 



the perturbed stellar surface. This can be achieved by imposing the continuity of the Lagrangian perturbation to the gravitational 
potential and its gradient at f = 1 : 

r nt + £ ■ V^ nt = *P ext + £ • W|f , W int + £ • V (V% nt ) = V^ 6 " 1 + £ ■ V (V^ xt ) , (12) 

where the superscripts "int" and "ext" correspond to "just below" and "just above" the stellar surface, respectively. Given that the 
gradient of the equilibrium gravitational potential is continuous, the first condition simplifies to 

xpint _ ,pext (B) 

The second condition can be simplified, using Poisson's equation, applied to the equilibrium gravitational potential: 

drT* = d^ n - 4 *? P °f, l V , d e T nt = deV ex \ <9/F int = <yF ext . (14) 
r 2 + ri 
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Here, ^ is the radial component of the Lagrangian displacement, when decomposed over the alternat e basis (a^, a g , a^) i ntroduced 
Reese et al. (2006 see also Eq. (A. 3 i). The first equation is different from what was applied in Reese et aT] ( 2009a i, since it 



takes into account the contribution from a non-zero surface density. However, given its low value, the resultant difference is quite 
negligible. The latter two equations are implied by Eq. ( fT3"] l, making it unnecessary to impose them. 

Finally, the usual mechanical boundary condition on the Lagrangian pressure perturbation, £p = 0, has been replaced by a 
slightly different condition: 



E ( ■ V 



(15) 



where E*> is a vector perpendicular to the surface (see Eq. (A. 2 1). The above, modified mechanical condition corresponds to setting 
the vertical gradient of 6p/Po, rather than its value, to zero. We note that Pesnell (19901 and |Dupret et al. (2002i applied a similar 
condition for spherically symmetric stars. In order to avoid having a boundary condition with a radial derivative in it, we calculate 
the dot product between Euler's equation and and cancel out the term corresponding to the vertical gradient of Sp/P{). This leads 



to a complicated expression which is given in spheroidal coordinates in App. A.2 



The main purpose in using Eq. ( p"5j ) is to obtain a non-zero value for 6T/Tq, useful for visibility calculations as described 
below. Indeed, when combined with a non-zero surface pressure and the adiabatic relation, the simpler condition, 6p = 0, leads to 
ST/Tq = 0, as illustrated in Fig.[T](left panels). One may then wonder if Eq. ( 15 i has an important effect on the frequencies and on 



the displacement at the surface. Numerically, it turns out the frequencies vary little when using either boundary condition, at least 
in the present study. The middle panels of Fig.[T]also show that the displacement is hardly affected. 
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Fig. 1. Lagrangian temperature variations and Lagrangian displacement at the stellar surface for two modes, using either dp = or 
Eq. ( fl3j ) as an external boundary condition. The figures to the right display the meridional cross-section of the Eulerian pressure 
perturbation of the two modes, divided by the square-root of the equilibrium density. 



2.6. Numerical method 

This set of equations and boundary conditions is projected onto the spherical harmonic basis and discretised in the radial direction 
before being solved using the code TOP ( |Reese et al. . 2006 2009a). Besides the improved way of treating the equilibrium model, 
the present calculations also benefit from a new form of finite differences. This form achieves 4th order accuracy for 1st order 
derivatives in spite of using windows with 4 rather 5 grid points. More importantly, this approach is robust to problems like mesh 
drift and spurious solutions. 



3. Mode visibilities 

At this point, we will switch to using the spherical vector basis (e r ,ee,e$) in which the polar or z-axis is lined up with the star's 
rotation axis. Furthermore, in order to make the equations more compact, we will prefer the notation "r" to "R s " when designating 
the stellar surface, although it should be understood as r{( = 1,9) in what follows. Similar implicit arguments also apply to other 
geometric terms such as r^ and rg. We introduce the unit vector e b s . which points from the star to the observer. Furthermore, we 
will assume that the vector e t, s . lies in the meridional plane <p = 0. Let i be the inclination angle, i.e. the angle between e, and e b s ., 
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where e, is lined up with the rotation axi^J An explicit expression for e b s . in terms of the usual spherical basis is: 
£obs. = (sin i sin cos cp + cos i cos 0)e r + (sin i cos cos <p - cos i sin 9)eg - sin i sin <pe,p 
The radiated energy, received by an observing instrument from a non-pulsating stars, is: 



E = r-^ ff /(y",geff,r eff )e obs . • dS 

M JJvis.Surf. 



(16) 



(17) 



where d is the distance to the star, /j the cosine of the angle between the outward normal to the surface and e b s ., g e s and T e ff the 
effective gravity and temperature, and /(ju, geff , 7W) the specific radiation intensity, multiplied by the instrument's and/or filter's 
transmission curve, and integrated over the wavelength spectrum,. As can be seen, the integral is carried out over the visible surface. 
In a pulsating star, this quantity is perturbed as follows: 



AE(t) = 



1 



27Td 2 



ff ' 

JJA(Vis.Surf.) 



Kp, ges , T eS )e ohs . dS + 



ff 

JJvis.Surf. 



Slip, g e s, T eS , t)e ohs . ■ dS + 



ff 



%,geff.7W>obs. -S(dS) 



Vis.Surf. 



(18) 

where (5 denotes the Lagrangian perturbation, 9v{. . . } the real part, and A(Vis.Surf.) the amount by which the visible surface is 
modified due the modifications of the surface normal induced by the oscillatory motions. Furthermore, we assume a complex form 
for the eigenfunctions, hence the reason for taking the real part of the above expression. 

[1977). 



The first term in Eq. ( 18 1 is proportional to the square of the displacement and therefore neglected (e.g. Dziembowski 



The Lagrangian perturbation to the specific intensity which intervenes in the second term may be developed as follows: 



81 = 1 



d\nl 5T e fi 
<91nr eff T eS 



din/ r5g eff \ | dl 

Slngeff geS I d/J 



(19) 



The partial derivatives of / are obtained from model atmospheres and will be dealt with in Sect. 3.4 We note that the above 
expression is not exact as it neglects the slight Doppler shifts caused by rotation and the oscillations. Such shifts modify the position 
of the emerging flux with respect to the instrument's and/or filter's transmission curve thereby modifying /, but are expected to play 
a negligible role compared to other effects. We are therefore left with a number of geometrical terms to calculate as well as the 
Lagrangian perturbation to the effective temperature and gravity. 



3. 1 . Geometrical terms 

In what follows, we will use the following expression for the displacement: 

£ = , 9, <f>, t)e r (0, 4>) + 6, 4>, t)e 9 (0, 4>) + 0, <f>, t)e, 



(20) 



where (e r ,ee,e^) is the usual vector basis in spherical coordinates. Note: we have used subscripts rather than superscripts, for the 
letters r, 0, and <p, to distinguish these components from those given in Eq. (A.4i. In what follows, we will be using the following 
relations: 

d e e r = e e , d e e g = -e r , d g e$ = 0, 

d$e r = sinf?e0, d$eg = cos6e$, d^e^ = - sin 0e r - cos 0eg 



(21) 



We start by calculating a surface element on a rotating, non-pulsating star. This is given by the following expression 

dS = (d s r x cfyrj d0d(/> = (r 2 sin 0e r - rr 6 sin Gee) d0d<f>. 

We then calculate the Lagrangian perturbation to a surface element: 

6(dS) = (dg^xd^,r + dgrxd^)d0dcf> 

= \^2r sin 6% r + r cos 0^ e + r sin 0dg%g + rd^A e r + (-r e sin 0% r - r sin 0dg^ r + (r sin - r g cos 6)£g - r s , 
+ y-rd^r + rgd^g + (r sin 9 - r e cos 0)fy) e } d6»d0. 



(22) 



eg 



(23) 



The approach used to obtain Eqs. ( p2| and (23 1 is essentially the same as that of Buta & Smith ( 1979} and Townsend ( 1997| l, but the 
effects of horizontal Lagrangian displacements are also included in the latter equation. We then calculate [i: 



H = <?obs. • « 



[r(sin i sin 9 cos + cos i cos 0) - ^(sin i cos cos <p - cos i sin 0)] 



V 2 + 4 



1/2 



(24) 



As opposed to a non-rotating star, the inclination cannot be arbitrarily set to to simplify the calculations. 
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where n 



AS 

HdS]| 



. The Lagrangian perturbation to yu is calculated as follows: 



f 6dS dS „ 1 f 6dS ( 6dS \ 

6, = e obs , ■ 5n = e obs . ■ |— - pjp*lldS||] = • [ m - [n ■ m jn 

[r e (sin i sin 9 cos <p + cos i cos 9) + r(sin i cos 9 cos cf> - cos z sin 0)] • [r fl ^ r - rdg% r + r£ fl + rgdg^g] 



V 2 + 3] 

sin i sin -rd^ r + rgd^g + (r sin - r e cos 
sinfl 



3/2 



(25) 



r[r 2 + r fl 2 ] I/2 

where we have used the relation £||dS|| = dS ■ £dS/||dS|| = n ■ 6dS. In the spherical limit, the above expressions become 
dS = (r 2 sin<9e r )d0d0, 

6(dS) = {(2r sin 0£ r + r cos 0ij g + r sin Qdoije + rd^A e r + (— r sin 0dg£ r + r sin 9%g) eg + (-rd^r + r sin 0%$) d#d</>, 
yu = sin / sin 9 cos ^ + cos i cos 0, 

[sin ; cos 9 cos - cos i sin 0] • \-dg( r + £g] sin ; sin (f> —d^r + sin 0& 

(5yU = • , 

r sinfl r 

in full agreement with the expressions previously obtained by |Heynderickx et al.| (fl994) for a non-rotating star with i = 0, provided 
one uses the relations: 

% r = 6r, & = r60, ^ = r sin 06<f>. (30) 



(26) 

(27) 
(28) 

(29) 



3.2. Lagrangian perturbation to the effective gravity 

Although the pulsation equations are given for a cylindrical rotation profile, we will relax this assumption in this section and allow 
for a general 2D rotation profile when establishing an expression for the Lagrangian perturbation to the effective gravity. Before 
giving the Lagrangian perturbation to the effective gravity, it is useful to recall various expressions for the unperturbed surface 
effective gravity: 



geff 



= -V¥ n + sQ. 2 e< = 



VP d( p o 



E? = - 



geffl. 



(31) 



Po Pi) 

The minus sign in the last equation comes from the fact that gravity is pointed inward, g e s being the norm of g e ff- Furthermore, 
contrarily to what happens in the stellar interior, the 9 derivative of Po vanishes at the surface since the surface is in pressure 
equilibrium. Any of the above expressions may be used to evaluate the surface effective gravity, although -VFo + sQ?e s provides 
the most accurate numerical results. It is also worth noting that the above expression neglects any contributions from meridional 
circulation and viscous forces to the equilibrium model {e.g. Rieutord & Espinosa Lara 2009 1. 
The relative Lagrangian perturbation to the effective gravity is then given by: 



Sgeff <%effll n ■ 6g eS 



\\g. 



effll 



geft 



(32) 



where we have used the simplification <5||g e ffll = ges • SgeslWge&W = ~ n • ^geff- The quantity 6g e f[ represents the vectorial Lagrangian 
perturbation to the effective gravity. It includes the Lagrangian perturbation to the gradient of the gravitational potential, and the 
acceleration of a particle tied to the surface, resulting from the oscillatory motions: 

S geB = -VT - £ • V (Wo) + (co + mfl) 2 f - 2i(co + m£l)Sl X £ - SI X (SI X £) . (33) 

Rather than working with the above expression, it is more useful to introduce the equilibrium effective gravity by adding and 
subtracting £ • V (sfi 2 e s ): 



S geif = - V¥ + £ ■ Vg eff + (co + mO) 2 £ - 2i(co + m£l)Sl x^-Slx(Slx^)-^-V (sQ 2 e s ) . 



(34) 



After some simplifications based on Poisson's equation and a lengthy derivation described in App. [B] one obtains the following 
explicit expression for a general rotation profile: 



(r 2 + rjjf' 2 d ( W 



S'eff 



rede* 



rr ( 



r r z + r: 



sin 9 (r sin 9 - rg cos 9) 



d ( {p}) + sin6» cos 9dg (il 2 ) 



, (rgdgg* (2r -r e cot 9) (r 2 + rj) + rr 2 e -r 2 r ee | ( rd e g eff rr e (~2r 2 - 3r 2 + rr g g) + (r 2 + r 2 ) rj cot 9 ) 
+ s> l ~5 T ; ^ 8eff } + & i o ; -3 geB I 



„2 i „2 



-(co + mQ) 



2 r^ r -rg^ e 



(r 2 + rl) 



1/2 



(r^r 2 ^ 
+ l -2i(u> + mQ.)n^ + 



r£ r - rg£g 



i0d ( (Q 2 ) + & sin 9dg (Q 2 ) 



(r^rl)' 
r sin 9 - rg cos ( 



(ri + rj) 



1/2 



(35) 
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It is worth noting that in the above expression, all of the terms involving a ( derivative of one of the perturbed quantities are divided 
by rf . The remaining terms involve no £ derivatives whatsoever. This characteristic is what one would expect for a quantity which 
is independent of the mapping. 

If a cylindrical rotation profile is used instead, the terms in square brackets simplify to the following expressions, respectively: 

sd s (fi 2 ) , (sin 0ft + cos 0ft) sd s (fi 2 ) . (36) 
In the non-rotating limit, Eq. ( [35| reduces to: 

SgtB = d r V + AnGp^ r - ^g eff - « 2 ft (37) 
r 

in full agreement with Dupret et al. (2002). 

Given the complexity of Eq. (|35|l, it is interesting to see if it can be approximated by a simpler expression. We consider the 
following approximation: 

6g™™- = -(u + man ■ n = -(*> + mQ) 2 (38) 

r + n) 

Figure [2] compares 6g e g and 6g^ rox ' for three different modes - two gravito-inertial modes and one p-mode. As can be seen in this 
figure, 5g^ mx ' is a very good approximation for p- and g-modes with a sufficiently high frequency. It is only for low frequency 
g-modes that the difference between the two expressions becomes non-negligible. 



cj/Q k = 2.0005 "/^x= 5.0400 u/Q K = 12,0905 




SO 40 00 60 SO 40 60 B0 30 40 60 SO 

Fig. 2. Comparison between 6g e ff and 6g^ lox ' for three modes. The upper row shows the meridional cross-section of the three 
modes, whereas the lower row compares the two calculations for the perturbation to the effective gravity at the stellar surface, 
between one of the poles (0 = 0°) and the equator (0 = 90°). For the middle and right mode, the two curves nearly overlap, making 
it hard to distinguish them. 



3.3. Lagrangian perturbation to the temperature 

We then turn our attention to the relative Lagrangian perturbations of the effective temperature, ST^g/T^. Given that the pulsation 
calculations are done in the adiabatic approximation, this quantity is not readily available. We will therefore content ourselves with 
the approximation 5T e slT s s m 6T/Tq. We note that this approximation is different from the approximation 6T e ff * ST given that 
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r e ff is latitude-dependant whereas To is not in barotropic stellar models. As was pointed out in Sect. |2.5| a modified version of 
the mechanical boundary condition was needed in order to obtain a non-zero value for 6T at the surface. One may then wonder if 
other approximations may be more suitable, such as applying the boundary condition dp = high up in the stellar atmosphere and 
extracting 6T at a deeper optical depth, or working with the Eulerian temperature fluctuation instead. The main difficulty with the 
first option is that the stellar models used in the present study do not include an atmospher^] Hence, it is not clear at what depth 
(which is likely to be latitude-dependant), one should extract the Lagrangian temperature variations. Using the Eulerian temperature 
variations also has problems of its own. Indeed, at the surface, the Eulerian temperature variations are dominated by the the advection 
term, thereby making them an order of magnitude larger than the Lagrangian variations, and with the opposite sign for the most 
part. Such variations do not seem the most appropriate from a physical point of view, since the intensity variations are more likely 
to result from local temperature variations physically present in the fluid, i.e. Lagrangian variations. Hence, given the limitations of 
our models and pulsation calculations, using the Lagrangian temperature variations along with the modified mechanical boundary 



condition seems like the best choice. Nonetheless, according to Dupret et al.|l|2002) and|Dupret et al.|([2003|l, adiabatic calculations 



of 6T/T0 are not reliable in the superficial layers and consequently lead to a poor approximation of 6T e g/T e ff. In order to obtain 
6T e ff/T e s accurately, one would need to do a full non-adiabatic calculations including a detailed treatment of the atmosphere such 
as what is done in Dup ret et al.|(|2002|l, usi ng 2D rotating models in thermal equilibrium, such as what is being developed in the 
ESTER project (Rieutord & Espinosa Lara, 2009, Espinosa Lara, 2010). 

The quantity 6T/Tq can be deduced very simply from 6p/Po via the adiabatic relation: 

ST_T 2 -16p m 



To r 2 p 

' lh] r > -.id 



where T2 is the second adiabatic exponent, given by the expression = y dl " ' 



3.4. Specific intensities 



Contrarily to the non-rotating case, the effective temperature and gravity of the equilibrium model depend on the latitude. 
Consequently, it is necessary to find the appropriate set of specific intensities and their derivatives for each latitude. We there- 
fore calculated a grid of K urucz atmosph eres with solar composition which spans the relevant temperature a nd gravity ranges. W e 
used the ATLAS 9 codtf] ( jKurucz] |1993|) in a modified form so as to include the convective prescription of |Canuto et al.| ( |1996| l, 



known as CGM (for more details see |Heiter et al.|2002| and [Barban et al.|2003| l. Due to convergence problems, some of the grid 
points were missing and had to be interpolated from neighbouring points. The final grid is illustrated in Fig. [3] (upper left panel) 
as is log I0 (geff) as a function of T e ff for a set of 2 M models with rotation rates ranging from 0% to 80% of the critical break-up 
rotation rate. 



For each grid point, we then approximated the p dependence through a Claret type law ( Claretj [2000 ) using a least-squares fit 



This yielded a set of five coefficients (I(p = 1), /(l)fli, /(l)fl2, ^(1)«3> /(1)«4) from which it is possible to calculate / and dl/dp 
for any p.. Using this approach rather than dealing directly with the original values for / yields better numerical results, both for 
interpolating / and especially for taking its p derivative. 

We then successively applied cubic spline interpolations, first as a function of log 10 (g e ff)> then as a function of T e s. This yielded 
b-spline coefficients from which it is possible to deduce the Claret coefficients and their derivatives for any value of log 10 (g e g) and 
r e ff within the relevant range. From this representation of /, we then calculated the Claret coefficients for /, dl/dg e g and dl/dT^ 
as a function of latitude for each stellar model. The remaining three panels of Fig. [3] display these quantities as a function of p 
and colatitude, 8, for a 2 M model at 80 % of the critical rotation rate. As is illustrated in these panels, both limb and gravity 
darkening are taken into account. The latitude-dependant Claret coefficients were then projected onto the spherical harmonic basis, 
truncated at 40 terms. This last step allows us to easily interpolate these coefficients onto a denser latitudinal grid during the visibility 
calculations. 

3.5. Visibility integrals 

Once the above quantities are calculated over the stellar surface, the visibilities may be evaluated numerically. However, rather than 



Eq. ( 18 1 into the following schematic form: 



using the time dependant form given in Eq. ( 18 I, it is more useful to extract the amplitude and phase of AE(t). To do so, we cast 



AE(t) = %\ if [A(6,4>) + iB(9,(/>)]e i "" t ' +i0J 'ded4) 

IJJvis.Surf. 



= % { \ \ [(A cos(m(f>) - B sm(m(f>)) + i (A sm(m(p) + B cos(m0))] e' w 'ddd(f> 



JJvis.Surf. 

cos(wf) I I (A cos(m<p) - B sin(m0)) d6d<f> - sin(«f ) J I (A sin(m^) + B cos(ot0)) dOdcp (40) 

JJvis.Surf. JJvis.Surf. 



2 The Kurucz atmospheres described in the following section are only used to calculate the emergent intensities as a function of r cff and g eff , 
and are not "joined" to the present models. 

3 See http : //kurucz . harvard . edu 
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Fig. 3. (Upper left panel) Effective temperatures and gravities of the Kurucz model atmospheres (symbols) and the 2 M e models 
(continuous lines; the lower longer lines correspond to more rapid rotation). The pluses (+) represent initial Kurucz atmospheres 
whereas the diamonds (o) represent missing points which were interpolated from the neighbouring atmosphere models. (Remaining 
panels) Specific intensities and its partial derivatives for a 2 M , 0.8 £2 C model, as a function of colatitude and /j. 



where A(6, <p) and B(8, (p) are real. This leads to: 

AE(t) =Ccos(wf + ip) (41) 



where C = V^t 2 + & is the amplitude and if — arctan(S/yi) the phase. Based on symmetry considerations with respect to the 
variable <p, the term & cancels out in the adiabatic case, thereby leading to C — \ and ip — or n. If, however, one tries to simulate 
non-adiabatic effects by introducing a phase-lag in the temperature variations, then the term S no longer cancels out. 

The integrals J{ and T> are evaluated numerically. The surface is discretised in the 6 direction, using a Gauss-Legendre collocation 
grid with typically Ng = 300 points, and in the <p direction using a uniform grid with a typical resolution of N$ = 720. The integration 
weights in the 9 direction are deduced from a Gauss quadrature, whereas those in the <p direction are uniform. Figure [4] shows the 
differences which result from increasing either Ng (left panel) or N$ (middle panel) in a 2 M Q stellar model, rotating at 60 % of the 
critical rotation rate. As can be seen, these differences are several orders of magnitude smaller than the visibilities themselves. 

At each point on the surface, the following condition is evaluated to determine whether or not the surface element is facing the 
observer and should be included in the integral: 

fi > (42) 



where /i is given in Eq. (24i. We note that this condition is only valid if no obstacles are present between the surface element and 
the observer, such as what could happen, for instance, in the more distorted SCF models with concavities at the poles (due to highly 
differential rotation) and an orientation which is not too far from equator-on. For such configurations, a more general condition 
based on a z-buffer or ray-tracing approach would be required. 

One possible concern is that using a Gauss quadrature implicitly assumes a spectral decomposition of the function which is 
being integrated, which could lead to a Gibbs phenomena near the cutoff between the visible and hidden side of the star. However, 
given that the Gibbs phenomena is oscillatory, one can expect its integrated contribution do be small. Nonetheless, we carried out a 
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test on the same model, using both Gauss quadrature and a trapezoidal integration. As can be seen in the right panel of Fig. [4] the 
two approaches yielded very similar results, showing that the Gibbs phenomena does not play an important role. 



8 resolution <fi resolution Int. method 




6 l\ 10 1 2 H 16 6 8 10 12 14 16 6 8 10 12 14 16 

7:rrc:ur::ic:y, Frequency, cj/Q k Frequency, cj/f) K 

Fig. 4. Visibility calculations and differences for a set of 3197 modes, obtained using two different resolutions, Ng = 300 and 
2000 (left panel), two different <p resolutions, N$ = 720 and 4000 (middle panel), and two different integration methods in the theta 
direction (right panel). These modes were calculated in a 2.0 M G model at 0.6 Q.c, with an inclination angle i = 10°. If i = 80° is 
used instead, the differences are reduced. 



3.6. Comparison with 1D case for non-rotating stars 



It is also important to check that in the non-rotating case, visibility calculations based on the above 2D integrations agree with the 
results obtained with a simpler ID integration, as is commonly used in other works (e.g. Heynderickx et al. 1994| l. In order to 
obtain the ID integrals, one starts with a given mode, which in this case will be proportional to a particular spherical harmonic, 
Yf(6, <p), re-expresses this spherical harmonic in terms of a new coordinate system with the z'-axis aligned with the line of sight via 



the following decomposition (Edmonds |1960| l 



17(0,0)= ^ d e mm ,(i)Yf(e',n 



(43) 



and calculates the integrated contribution from each of the Yf (ff , (/>'). In this new coordinate system, the visible surface corresponds 
to a hemisphere. Consequently, all of the non-axisymmetric terms m' + cancel out and only the term d e m0 (i)Y®(8' ,$') remains. 
This leads to the following expression: 



AE(t) 



d 1 



T e ff 



f 

Jo 



ai 



<91nr eff 



Yfc'dfx' + 



eff r 

ff Jo 



31 



ges Jo dlngeff 



y&w + -(2+m-t) 



f IYfa'dfi' 
Jo 



(44) 



where 



4o(0 



47T 

27TT 



Yf (6 = i, 



0), 



cos ( 



ST, 



elf 



and 



?eff = (<5|eff) IT 



and s — +1 depending on the sign convention used in defining the spherical harmonics, and the values of i and m. Figure [5] shows 
the differences between the two approaches. As can be seen, the results are very similar, thereby validating the 2D integrations. 



3.7. Comparison with Lignieres et al. (2006) and Lignieres & Georgeot (2009) 

Before going on to describe the results, we carry out one last test in which we compare the 2D integration method described here 
with the harmonic projection method used in Lignieres et al. (2006} and |Lignieres & Georgeot ( 2009| l. In order to carry out the 
comparison, we calculate simplified disk-integration factors, D(i), by only integrating the temperature fluctuations over the visible 
surface. We avoid using the intensity fluctuations deduced from the Kurucz atmospheres, as this would introduce the effects of limb 
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darkening and a <p dependence of a different form than e' m ^, which is not currently implemented in the method by Lignieres & 
Georgeot (2009). Hence, the disk-integration factors are given by: 



D(i) cos(wf + if/) 



1 



I nRL (6T) JJvis.Surf. 



6Te n 



as 



where 



(ST) 



411 



\6T\ 2 dS 



1/2 



(45) 



and where S is the total surface and ift a suitably chosen phase. The disk-integration factors are normalised in the same way as 
is done in Lignieres et al. ((2006) and Lignieres & Georgeot (2009), which leads to D(Q) = 1 if ST = 1. Figure [6] shows the 
differences between the two approaches for a model rotating at 0.8 Qc- As can be seen, the differences are quite small relative to 
the visibilities themselves. At this point, it is helpful to remember that the curve which delimits the visible surface, is approximated 
as the intersection of a plane and the stellar surface in the harmonic projection method. One can expect such an approximation to 
become less good at high rotation rates due to the strong centrifugal deformation, but as shown in Fig.[6]the effect remains small, in 
full agreement with the results presented in Lignieres & Georgeot ( 2009), based on a Roche model. 
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Fig. 5. Visibility calculations and differences for a set of 
1223 modes in a non-rotating model, using ID and 2D in- 
tegrations. The latitudinal resolution is Ng = 300 for both 
approaches, and N$ = 720 points are used in the azimuthal 
direction for the 2D integrations. The inclination angle is 
i = 10°. Note: because of mode degeneracy, there are much 
fewer frequencies. 



Fig. 6. Simplified disk-integration factors and differences for 
a set of 1507 modes in a model rotating at 0.8 Qc, using 2D 
integrations and the harmonic projection method described 
in Lignieres et al. (2006 ) and Lignieres & Georgeot (2009). 
The resolutions are Ng = 300 for both approaches, and 
N,f, = 720 for the 2D integrations. The inclination angle is 
i = 10°. Note: due to symmetries in the simplified visibility 
calculations, the visibilities of some modes cancel out and 
were therefore not included in this plot (see text for details). 



Given the simplified form given in Eq. (45 I, it turns out that modes with |m] > 2 and such that (— l) m+ 8T(n - 0) — 8T{6), i.e. 
even modes when m is odd and vice-versa, have zero disk-integration factors, regardless of orientation. Consequently, these lead to 
numerical results around 10~ 16 or less, for both methods, and were therefore not represented in Fig. [6] Appendix [c] explains why 
these disk-integration factors cancel out. 



4. Overall visibilities 

In what follows we will first focus on mode visibilities in the CoRoT photometric before dealing with amplitude ratios between 
different bands. 



4.1. General results 

Fig.|7]shows the visibilities of modes within acoustic frequency ranges for a set of 2M Q stellar models at 9 different rotation rates 
(1 for each row). The columns correspond to 4 different inclinations. Given that the intrinsic mode amplitudes cannot currently be 
predicted, as this would require a full non-linear development, we normalised the modes so that the maximal displacement within 
the entire star, ||£|| max , multiplied by the square of the co-rotating frequency, (o> + mQ.) 2 , is constant. This tends to favour acoustic 
modes, for which the maximal amplitude is near the surface. Multiplying by the square of the frequency yields a near constant value 
for the visibilities of modes with similar I and m in the acoustic asymptotic regime, for the non-rotating model (upper row). The 
differ ent colours indicate the I value w hich woul d carry on fro m Q 
2006: 



In 



Lignieres & Georgeot (2008) (see also 



Reese 



through a mode-following procedure (e.g. Lignieres et al. 



2008 ), a different set of quantum numbers (n, 



m) is introduced. This set 

of quantum numbers is better adapted to island modes, the rapidly rotating equivalent of modes with low I - \m\ values. The number 
h corresponds the number of nodes along the underlying ray path associated with these modes, whereas £ is the number of nodes 
transverse to this ray path. In the last row of Fig. [7] the visibilities of the 0.8 £2c model are repeated, but using a colour scheme 
which indicates the value of t. 
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Fig. 7. (Colour online) Mode visibilities for 2 M stellar models at 9 different rotation rates, and 4 different inclinations (i = 0° 
corresponds to a pole-on configuration). The colours correspond to the I value carried on from Q = 0, except for the last row where 
they correspond to I. 
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A first trend which appears at rapid rotation rates is that the equator-on configurations tend to have larger amplitudes. This is 
simply due to the fact that we normalise the visibilities with respect to the stellar luminosity in the observer's direction. Given that 
the star is less luminous equator-on, this leads to larger relative amplitudes. If instead we normalise by the pole-on luminosity, such 
as what is done in Lign ieres et al.| ( |2006| l, the amplitudes show a much smaller increase, due to the equatorial focusing of island 
modes, as will be discussed in the following section. 

A second effect which can be seen in these panels is that the pulsation spectra become more and more complex as Q increases. 
This is caused by several factors. First of all, as opposed to the non-rotating case, modes with the same n and I values but different 
m values have different frequencies. At small (uniform) rotation rates, this leads to frequency multiplets which are evenly spaced 
and which span a small frequency range. As the rotation rate increases, these multiplets become uneven and they start to overlap, 
thereby leading to a complex spectrum. Secondly, modes with intermediate ( - \m\ values become chaotic modes at rapid rotation 
rates. This causes their node spacing to become uneven and is likely to make cancellation effects less efficient than in the non- 
rotating case (Lignieres & Georgeot, 2009). Accordingly, the number of modes which are visible above a given threshold increases. 
One may wonder if this can explain the irregular nature of pulsation spectra in 5 Scuti stars. Up to some extent yes, but not entirely. 
In both cases, the frequencies are distributed so that no simple regularities stand out. However, observations show a somewhat 
larger dispersion in amplitudes, with only a few high-amplitude modes. This suggests that non-linear effects lead to intrinsic mode 
amplitudes which cannot be described in a simple way as is done here with the ad-hoc normalisation, thus making it difficult to 
identify modes from visibilities alone. 

Another important effect is the relative simplicity of the frequency spectra for the pole-on orientation, and to a lesser extent for 
i = 30°, even in the most rapidly rotating models. Indeed, non-axisymmetric modes cancel out in the pole-on configuration, and the 
I = modes, i.e. the rotating counterparts to t — and 1 modes, stand out. This means that such an orientation facilitates interpreting 
the oscillation spectrum of rapidly rotating stars. Of course, it is also more difficult to identify pole-on stars as rapid rotators, since 
spectroscopic observations will only reveal narrow absorption lines due to the unknown sin i factor (although the overall shape of 
the spectrum will differ substantially from that of a black-body at a single temperature). Interferometric studies have, however, been 



able to confirm that Vega is a rapidly rotating star, seen nearly pole-on (i = 4.54°, Peterson et al. 2006 1. This recently motivated a 




search for pulsation modes in this star. Bohm et al. (2012} found some stellar variability, but were unable to confirm the presence of 
stellar oscillations, thereby showing the need for further observations with larger instruments. It is also necessary to search for other 
pole-on rapid rotators to see if they exhibit stellar pulsations. 

In the equator-on configurations, the visibilities of the chaotic modes increases relatively to those of the island modes, thereby 
making it harder to distinguish the two. This result partially agrees with the results presented in Lignier es & Georgeot| 
although we do note that here, some of the island modes still tend to remain more visible whereas in Ligni eres & Georgeot| ' 
the two had comparable amplitudes. To understand where this difference comes from, we tried various manipulations. In 
we show alternate ways of calculating or normalising the visibilities for the model rotating at Q. — 0.8 Qc- The top row shows the 
visibilities of all of the modes, whereas the bottom row corresponds to axisymmetric modes (m = 0) only. The first column uses 
the original normalisation, so that the upper left panel of Fig. [8] is identical to the lower right panel of Fig. [7] The second column 
shows shows what happens if the modes are normalised by (dT^) rather than « 2 ||£|| m ax- The third column shows what happens 
if we use the approximation 5T 6T e g instead of 6T/T ST^/T^. Finally, the last column shows the disk-integration factors 
as given in Eq. (45 1. For the sake of legibility, we divided the visibilities in each panel by the maximal value, so that the results 
are between and 1 . As can be seen, the most visible modes still tend to be island modes regardless of how the visibilities are 
calculated and whether or not only axisymmetric modes are kept. Hence, another explanation is needed for the difference between 
the present results and those of Lignieres & Georgeot (2009). One difference between the two studies is the boundary condition on 



the pressure perturbation. Indeed, in the present study, we set the vertical gradient of dp/Po to zero at the surface (Eq. 15 I, in order 
to avoid having 6T — at the surface. Lignieres & Georgeot ( 2009 1 set dp = instead. This however does not lead to 6T — 
in their case, since they are dealing with a polytropic model, where the pulsation equations become singular at the surface. This 
difference in boundary condition may have an important impact on the temperature fluctuations at the surface which dominate the 
mode visibilities presented here, although it is hard to test since applying the condition 5p — with SCF models would only lead to 
ST = 0. One possibility would be to apply the condition Eq. ( 15 1 to polytropic models. 

Finally, Fig. [9] shows the separate contributions coming from the g e e, T e ff and dS variations. As can be seen in the figure, the 
effective temperature fluctuations are the most dominant effect. This is because effective gravity plays a small role in the specific 
intensity emitted by an atmosphere. The surface deformation plays a small role in acoustic modes because it is proportional to the 
displacement rather than the displacement times the square of the frequency. 



4.2. Equatorial concentration 

One of the consequences of the centrifugal force on acoustic modes is the equatorial focusing of island modes around the equator, 
due to the geometric distortion of the resonant cavity. Lignieres et al. ( 2006 ) pointed out the possibility that this could lead to higher 
observed pulsation amplitudes in stars observed in an equator-on configuration, thereby explaining the observations of |Suarez et al.| 
(2002), which show a correlation between mode amplitude and inclination. It is also worth noting that for gravito-inertial modes 
in the inertial regime, modes become trapped in an equatorial waveguide (e.g. Dint rans & Rieuto"rd| |2000| |Townsend] [2003b). 
Paradoxically, this usually causes these modes to be more visible in a pole-on configuration, due to the contribution from surface 
distortion (Townsend 2003a Daszyhska-Daszkiewicz et al. 2007) >. 

Figure 10 shows the visibilities of various axisymmetric (m = 0) acoustic modes as a function of inclination. The left panel 
corresponds to a non-rotating model, whereas the three middle panels corresponds to a model at 0.8 Oc- At such a rotation rate, all 
of the modes in these panels are island modes. The mode visibilities in the second panel are normalised by the pole-on equilibrium 
luminosity of the star (as is done in Lignieres et al. 2006). What is surprising in this figure is that although the symmetric modes 
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Fig. 8. (Colour online) Mode visibilities calculated or normalised in alternate ways for the model rotating at Q = 0.8£2c> as viewed 
from the equator. The upper shows the visibilities of all of the modes whereas the lower panel only includes axisymmetric (m = 0) 
modes. See text for details. 




Fig. 9. (Colour online) Separate contributions to the visibilities, coming from the geg, T e g and the dS variations. These are calculated 
for the model at Q = 0.8 Qc (which corresponds to the before last row of Fig. [7]). 



are more visible from the equator than from the poles, one would expect a higher contrast between the two, especially for the I — 
modes. Indeed, the surface temperature fluctuations, 5T/T, of the £ = modes, depicted in the rightmost detached panel, show a 
much stronger confinement towards the equatorial region. A careful investigation reveals that the reason for this lack of contrast 
between polar and equatorial visibilities stems from the geometrical shape of the star. Given the high curvature at the equator, 
the surface quickly starts pointing more towards a polar direction. To illustrate this point more clearly, the third panel shows the 
same visibilities assuming the star is spherical, but retaining the other effects of rotation such as gravity darkening. Clearly, a 
much stronger contrast appears between the poles and the equator. One can also quantify at what latitude a temperature fluctuation 
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contributes more to mode visibilities in an equatorial direction than in a polar one by calculating at what point the ratio, 

2 Cf , 2 ■§- [n (6o, <p) ■ e x , g eS (9 ) , T etf (#„)] e x ■ dS 

r «*>) = t f r . (46) 

exceeds one. In the above formula, only the effects of temperature fluctuations are retained as they usually are the dominant factor in 
mode visibilities. The denominator corresponds to the intensity fluctuation in the polar direction, e : , caused by a unit axisymmetric 
effective temperature variation at colatitude 6q. The numerator is the intensity fluctuation in the equatorial direction, e x , caused by 
the same unit effective temperature fluctuation, and it's mirror image in the southern hemisphere, at colatitude n - Oq, hence the 



factor 2. In the rotating model used in Fig. 10 the ratio R(8q) exceeds 1 from 73.1°, which is depicted by the vertical line in the 
rightmost panel. As can be seen in the figure, this line is quite close to the maximum of the (n, () = (9, 0) mode, thereby explaining 
the near equal pole-on and equator-on visibilities. Had the star been spherical, the position of this line would have been 57.0°. 



Non^rotatmg AE/E pole _ on AE/E observer 6T/T 

0.20{ + t t 1 I 




20 40 60 800 20 40 60 800 20 40 60 800 20 40 60 80 20 40 60 80 

Inclination, i (in degrees) 9 (in degrees) 



Fig. 10. Visibilities of various island modes as a function of the stellar inclination. The t and n values are indicated by the styles and 
shades of the lines. The rightmost detached panel shows the temperature fluctuations, ST/T, of the two I — modes in the rotating 
model. 



The fourth panel shows the mode visibilities if they are normalised by the equilibrium luminosity in the observer's direction (as 
is done in Fig. Once more, there is a strong contrast between the polar and equatorial visibilities. This is simply caused by the 
lower luminosity of the star equator-on, which results from gravity darkening. From an observational point of view, this last ratio 
is easier to obtain given that £p i e -on is not directly observable, and may be one of the dominant factors in the coiTelation between 



inclination and mode amplitude found by |Suarez et al. (20021 



4.3. Avoided crossings 



It is then interesting to look at the effects of avoided crossings on mode visibilities. Lignieres et al. (2006 1 pointed out the possibility 
that rotationally-induced avoided crossings may explain close frequency pairs observed in 6 Scuti stars ( jBreger & Bischof] [2002 



Breger & Pamyatnykh 2006 ). Indeed, the structure of the eigenfunctions of the two (or more) modes involved in an avoided crossing 



are a mixture of the eigenfunctions prior to the crossing. Hence, if a high-visibility mode interacts with a low-visibility mode, the 
two may be visible during the crossing when both modes have a mixed character, as was already shown in Daszyriska-Daszkiewicz 

[etaLl ( [2b02l>. 



Figure 11 (left panel) shows a set of avoided crossings between families of I — 1, 5 and 9 modes, the azimuthal order being 



m — 1 . As can be seen in the figure, rotation is an ideal way of generating avoided crossings due to its differential effect on modes of 



different degrees. As was found in Lignieres et al.|(|2006|l, we also find that modes of degree I tend to interact with modes of degree 



I + 4. The middle panel shows the associated mode visibilities. At 0.48 Qc, the modes labelled (a) and (b) are strongly interacting, 
thereby producing an intermediate visibility for both modes. We note that the frequency separation between the two is 2.24 //Hz, 
which is slightly larger than for observed close frequency pairs. The meridional cross-section of these modes is depicted in the right 
panel, both of which show mixed characteristics. 

In order for avoided crossings to explain close frequency pairs, a number of conditions need to be meet. On the one hand, the 
coupling between the modes needs to be sufficiently weak so that the frequency differences become small enough. In the above 
example, none of the frequencies get closer than 1 .4 //Hz which is slightly above what is considered to be a close frequency pair. 
On the other hand, the coupling needs to be sufficiently strong so as to allow the modes to have mixed characteristics over a non- 
negligible range in parameter space. Then, one would need to do a full statistical study of the frequency spectrum of visible modes 
and quantify how many close frequencies are due to avoided crossings and how many result from chance coincidence between 
uncoupled modes. Such a study is, however, beyond the scope of this paper. 
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Frequencies Visibilities 




Fig. 11. (Colour online) A set of avoided crossings between £ = 1, 5 and 9 modes with m = 1. The linestyles and colours indicate the 
I value at Q. — 0.3 Qc (note that an avoided crossing has already taken place between the uppermost t = 5 and 9 modes, hence the 
"swapped" identification compared to the other modes). The left panel shows the frequencies, the middle panel gives the visibilities 
in the CoRoT photometric band for i = 30° and the right panel shows the meridional cross-sections of two modes undergoing an 
avoided crossing, and labelled by "(a)" and "(b)" in the two left panels. 



5. Multi-colour visibilities 

One of the advantages of calculating amplitude ratios between different photometric bands is that the intrinsic amplitude of the 
mode factors out. As a result, it is possible to predict such ratios from linear theory alone. In this section, we will examine some 
of the effects of rotation on amplitude ratios. In particular, we will focus on a frequency multiplet, and on modes with the same I 
and m values. However, before discussing different results, it is useful to examine more carefully the question of normalising mode 
visibilities in multiple photometric bands. 



5.1. Normalisation 

In what follows, we will depart from the traditional approach which consists in normalising the visibilities with respect to a given 
band. The problem with such an approach is that sometimes, there are modes for which their visibilities nearly cancel out in the 
chosen band, thereby amplifying their amplitude ratios. This can then exaggerate differences between these ratios and those of other 
modes. We therefore propose another normalisation scheme which does not favour one band over the others, and which applies to a 
set of modes which we want to compare. 

We start by considering a set of N modes, numbered j — 1 . . . N. Let us denote Vj the visibility of mode j in the photometric 

band i. We will also use the notation V' to represent the vector (V/),=i...i. We want to normalise each mode so as to minimise, in 
some sense, the distance between their visibilities. To do so, we introduce a set of references visibilities Wj, which at this point are 
unknown, and normalise each mode so as to minimise its distance to W,. This amounts to minimising the following cost function: 

N b 

/= ZZ( A ; y /- w '') 2 ' (47) 

where b is the number of photometric bands and Aj the normalisation factors for each mode. In this expression, both the Aj and W, 
are unknowns. Furthermore, we impose the following constraint, so as to avoid the trivial solution Aj = 0, W, = 0: 



1=1 



2 - 1. (48) 



As will be shown in App. 



V = (Wi)i=i...b turns out to be the principal component of the normalised vectors (y J ^j ,_ ] ^, where 
Vj = V] I ^Ylk=\ iyii) ■ This result is fairly intuitive, since the vectors (t^) . j N retain the directional information (in b dimen- 
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sions) associated with each mode while discarding the amplitude (which is arbitrary anyway). Once the reference visibilities are 
determined, it is straightforward to find the normalisation factors: 

yb yJyy. b 

Aj = 1=1 ' ' , or equivalently, Aj = V V/W/. (49) 

xf , (vv) 

The normalised visibilities then become V J - AjV' — AjV 1 . This last equation has a simple geometrical interpretation: indeed, 
AjV J represents the projection of W onto the direction VK 

In the following sections, we plot normalised visibilities V J based on the set of modes relevant to each sub-panel of the various 
figures, and will call them "amplitude ratios" even if strictly speaking they are not amplitude ratios. In some cases, the reference 
visibilities, Wi, are also displayed. It is important to bear in mind that this normalisation depends on the set of modes, and furthermore 
on their inclination since this affects their visibilities in rotating stars. Hence, the mode-dependant normalisation factors will change 
from one panel to another, even for the same set of modes, if viewed at different inclinations. Arguably, some of the information 
based on inclination is lost in the process (the same would also be true if we had normalised the visibilities with respect to a given 
band). However, what is important here is the comparison of amplitude ratios between a set of modes for a given configuration. 
Furthermore, one does not have the luxury to modify the inclination of an observed star. Another concern is that by normalising 
the visibilities by a mode-depend scale factor, one is seemingly no longer comparing amplitude ratios which are the true invariants. 
However, such information is contained within the shape of the curve defined by the normalised visibilities, i.e. the ratios between 
the components (V/), = i.../,. Moreover, with the above normalisation, modes with similar amplitude ratios will have similar scale 
factors, Aj, and hence normalised visibilities which are similar, even more so than if one of the photometric bands had been used 
for the purposes of normalisation. 

5.2. A multiplet 

We first start by looking at the (n, €) = (6, 2) multiplet. Figure [12] shows amplitude ratios, based on the Geneva photometric system, 
at different rotation rates and inclinations. Only axisymmetric modes are displayed for the pole-on configurations (z = 0°) and even 
modes, i.e. modes which are symmetric with respect to the equator, are shown for the equator-on configurations (i = 90°), given that 
the other modes cancel out. As expected, the amplitude ratios do not depend on the inclination or azimuthal order in the non-rotating 



case. When the rotation rate increases, this is no longer true, as has been already pointed out by |Daszynska-Daszkiewicz et aL 



(2002 1 and Townsend (2003a). It is interesting to note that this departure from a constant behaviour is non-monotonic and therefore, 
somewhat difficult to predict. For instance, at (Q, i) = (0.4 Qc, 30°) the \m\ = 2 modes have higher ratios at lower wavelengths and 
lower ratios at higher wavelengths. At (Q, z) = (0.5 Qc, 30°) the opposite is true. The |m| = 1 modes have the opposite behaviour. The 
dependence on orientation seems to be somewhat simpler, although one can note how the mode order frequently changes between 
i = 30° and i - 60°. This implies that mode identification from amplitude ratios will not be as straightforward as in the non-rotating 
case. 

It is also interesting to note how the overall shape of most of the ratios change starting from Q = 0.6 Qc- At lower rotation 
rates, the concavity seems to be pointing up and the largest ratios are obtained at the lowest wavelength. In the most rapidly rotating 
models, the concavity is pointing downwards, and the ratios at the lowest wavelength are among the smallest. A particularly striking 
example of a change in behaviour is the drastic modification of the m = amplitude ratios between Q = 0.5 Qc and Q = 0.6 Qc, 
especially for the pole-on configurations (i = 0°). In order to gain a better understanding of what is causing this, we plot in Fig.fTT 
the meridional cross-sections of this mode for a selection of rotation rates. As can be seen, an avoided crossing takes place between 
Q = 0.5 Qc and Q = 0.6 Qc- This agrees with Daszynska-Daszkiew icz et aL| ( |2002) > who had also shown that avoided crossings can 
have an important effect on amplitude ratios and phase differences, and who furthermore showed how these can lead to erroneous 
identifications. A number of other avoided crossings involving the non-axisymmetric modes also take place around that rotation 
rate. Another effect which comes into play is the fact that the modes depart from the polar regions to focus around the equator 
starting from Q 0.6 Qc, as can be seen in Fig. [13 for the m = modes. 



One last interesting feature in Fig. 12 is the simi 



arity between prograde and retrograde modes with the same \m\ value. Townsend 
( |2003a| l had previously noted a similar phenomena for gravito-inertial modes, using the traditional approximation. He was, however, 
unable to provide an explanation for this behaviour given that the Coriolis force plays a dominant role in the modes he studied, 
thereby leading to important differences between prograde modes and their retrograde counterparts. The modes shown here are 
acoustic modes, in which the Coriolis force only plays a minor role. Hence, the structures of prograde modes and their retrograde 
counterparts are quite similar, thus explaining the similar amplitude ratios. 

5.3. Modes with the same I and m values 

Figure [14] shows the amplitude ratios for a set of modes with consecutive radial orders, n, and with {(, m) - (1, 1). In the non- 
rotating case, amplitude ratios for fixed (t, m) values are expected to be similar, though not identical, since the ratios between | r /r, 
<5Zeff/r e ff, and 5g e ff/g e s can vary slightly from one mode to the next. This is different from the purely geometrically disk-integration 



factors calculated in Lignieres et al. (2006 ) where strict equality is expected (and obtained) at Q = 0. What is interesting to note here 



is that for the most part, this similarity between the different ratios continues up to rapid rotation rates, as can be seen in Fig. 14 
The reason for this behaviour is straightforward. Modes with fixed (i, m) values produce island modes with the same {t, m) values 
These modes have a similar lateral structure perpendicular to the underlying ray path, apart from the width which decreases with 
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Fig. 12. Amplitude ratios of the (n, {) = (6, 2) multiplet, based on the Geneva photometric system. Each column corresponds to a 
different orientation and each row to a different rotation rate. The numbers on either end of each plot and connected by dotted lines 
indicate the azimuthal order, m. The continuous, dotted, and dashed linestyles correspond to m = 0, \m\ = 1 and \m\ = 2, respectively. 
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Fig. 13. Meridional cross sections of the (n, €, rri) = (6, 2, 0) mode at different rotation rates. At the lowest rotation rates, the modes' 
structure is dominated by the central parts. 



frequency, and only differ by the number of pseudo-radial nodes. Recently, Pasek et al. ( 2012 ) came up with asymptotic expressions 
for their structure, based on ray-dynamics. One may then wonder why we didn't directly select modes according to their ((, m) 
values. The problem with selecting modes that way is that it mixes together even and odd modes, which have different visibilities 
when the star is close to equator-on. Although the lateral structure is the same, the underlying ray path curves around the equator so 
that in odd modes, one of the pseudo-radial nodes actually corresponds to the equator. Hence, even and odd modes with the same 
{(, rri) values will have a similar surface structure, if one limits themselves to one of the hemispheres. Accordingly, one set of (t, m) 
values correspond to two different sets of {£, rri) values, one for even modes and the other for odd modes. 

It would, however, be incorrect to conclude that modes with the same t and m values systematically lead to similar amplitude 
ratios. In Fig. 15 we represent the amplitude ratios for modes with {£, m) - (3,0). As opposed to Fig. 14 this is one of the 
worst cases, where agreement is rather poor. Various effects can lead to differences in amplitude ratios. As mentioned above, 
avoided crossings can cause irregular behaviour. In some cases, a mode can be at the transition where temperature variations start to 
dominate over surface deformations. This can lead to a very different behaviour from the modes before or after the transition, and 
cause some of the visibilities to nearly cancel out. A good example of this is the n — 6 mode for i = 0° and D. = O.lfic- Finally, if 
the radial order is too low, a given mode may depart sufficiently from the asymptotic regime to produce noticeable differences on 
the visibilities. 



6. Conclusion 



In this paper, we derived a new set of equations to calculate mode visibilities in rapidly rotating stars. These equations take into 
account the centrifugal deformation of the star as well as gravity darkening, and treats the modes in a more realistic way by 
including the associated effective temperature and gravity variations, as well as the surface distortions. The modes are calculated in 
fully deformed 2D stellar models based on the SCF method (Jackson et al. 2005 MacGregor et al. 2007), using the 2D oscillation 
code TOP (Reese et al. 2006 2009a). A grid of Kurucz atmospheres was used to calculate realistic intensities at each point on the 
surface ( |Kurucz| |2005 1. As such, this represents an important step towards obtaining realistic multi-colour visibilities of acoustic 
modes in rapidly rotating stars, and an improvement over previous studies which approximated the effects of rotation on the model 
and the oscillations (e.g.iiDaszynska-Daszkiewicz et al. 2002 Townsend 2003a[ l or used simplified visibility calculations (Lignieres 
&Georgeot||2009] l. 



One of the important limitations in the present study is the use of the adiabatic approximation. Accordingly, the temperature 
variations are not reliable near the surface and, hence, provide a poor approximation of the effective temperature variations (Dupret 
et al. 20021. Including non-adiabatic effects is likely to have an important impact on mode visibilities since effective temperature 
variations often play a dominant role in mode visibilities, as was shown in Fig. [9] Nonetheless, one can hope that the present 
calculations will give insight into the qualitative behaviour of mode visibilities and amplitude ratios in rapidly rotating stars. In 
future studies, we plan to implement non-adiabatic effects, including a full treatment of the stellar atmosphere, by applying a similar 
approach to |Dupret et al.|(|2002|l. This will re quire the use of thermally relaxed rapidly rotating models, such as those produced by 
the ESTER code (Rieutord & Espinosa Lara 2009 Espinosa Lara 2010 1. 

The results presented here confirm a number of results previously established. For instance, in stars observed pole-on, some of the 
island modes stand out compared to other modes and form a regular frequency pattern (Lignieres & Georgeot 2009). Consequently, 
rapidly rotating pulsating stars seen pole-on may be promising asteroseismic targets. Chaotic modes are more visible than their non- 
rotating counterparts, probably as a result of irregular node placement, thereby complicating the frequency spectrum, especially for 
equator-on configurations (Lignieres & Georgeot |2009| ). Avoided crossings can substantially modify mode visibilities and cause 
modes which are not usually visible to exceed the detection threshold ( |Daszyriska^ Daszkiewicz et al. 2002 ). Lignieres et al. ( 2006 ) 
previously pointed out that these may explain observed close frequency pairs in 6 Scuti stars, although a full statistical study is 
required to test this possibility. Finally, amplitude ratios depend on both the inclination, i, and the azimuthal order, m, in rotating 
stars (Townsend 2003a Daszyriska-Daszkiewicz etal. 2007[ ). This, of course, makes mode identification in such stars more difficult. 

New results were also obtained. In particular, although rapid rotation causes island modes to focus around the equator, the 
corresponding visibilities show a much smaller contrast between pole-on and equator-on configurations (unless normalised by star's 
luminosity in the observer's direction). This is because the same geometrical distortion which causes mode focusing in the first 
place, also increases the proportion of the stellar surface which points in a more poleward direction. We also showed that acoustic 
modes with the same ((, \m\) values tend to have similar amplitude ratios, although this effect is not systematic. The similarity 
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Fig. 14. Amplitude ratios for sets of modes with consecutive radial orders and for which ((, m) = (1, 1). The reference visibilities 
used for normalisation, W,, are shown as a dashed line in each plot. In most cases, this line is covered up by the amplitude ratios, 
which are represented by continuous lines. The numbers on either end of the plots indicate the radial orders, n. 
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Fig. 15. Same as Fig. 14 except that (€, m) = (3,0). 
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between prograde and retrograde modes stems from the small influence of the Coriolis force on acoustic modes. Modes with the 



same {I, m) values have a similar surface structure, as expected from asymptotic ray theory (Pasek et al. 2012 1. 

In a forthcoming study, we plan to study the statistical properties of the frequency spectra by analysing their autocorrelation 
functions (Ligniere s et al.||2010| l and by looking at the cumulative distribution functions of the frequency separations. These will be 
compared with observations, in order to assess up to what extent the present theory is realistic and to see if it is possible to extract 
global quantities such as the large frequency separation or rotation rate from the observations. Using multi-colour photometry, we 
will also develop a new strategy for constraining mode identification and obtaining or confirming the values of global quantities. 
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Appendix A: Equations in spheroidal coordinates 

In this section, we derive explicit expressions for the pulsation equations and the mechanical boundary condition, based on the 
coordinate system described in Sect. |2.4| However, before giving these expressions, it is useful to recall a few definitions. The 
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natural covariant basis, denoted (Eg, Eg, E$), is defined via the relation E, = <9,r, where i stands for 8 or <f>, and r = re r : 

E ( = r ( e r , E g = r g e r + re g , E^ = rsin8e^. (A.l) 

Here (e r , eg, e$) is the usual spherical basis associated with the spherical coordinates (r 8, <p). The associated dual (contravariant) 

basis is defined such that E' ■ Et — 5'.: 

J j 



E? = 



E = 



r ( rr ( 



r sin 6* 



(A.2) 



The vector is perpendicular to surfaces of constant f value, including the stellar surface. As in Reese et al. (2006 1, we derive an 
alternate basis from (Eg, Eg, E$) as follows: 



r 2 rg 



r 2 rg 



dg - ~i—Eg, 

I- 



r 2 r^ sin 8 



(A3) 



In the spherical limit, the alternate basis converges to the spherical basis. The Lagrangian displacement is decomposed over the 
alternate basis as follows: 



t = $a c + faB + pa 4 



These components are related to the spherical components (see Eq. (20l) as follows: 



(A.4) 



(A.5) 



where superscripts are used with the alternate components, and subscripts with the spherical components. Based on the alternate 
components, the dot product £ • g e s becomes: 



% ■ geff = 



where we have introduced the following quantities: 
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A. 1. Pulsation equations 

We now give explicit expressions for the pulsation equations in spheroidal coordinates. The continuity equation is: 



0=* + -f 
Po r z r ( 



Euler's equation takes on the following form: 
= [co + mQ.] 2 
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Qf 2 sin0 . (rg sin 8 + r cos 8) „ 
if — 2i [o> + mil] £ 



po sin 8\Poj sin 8 sin 8 sin 8 
where s — r sinfl is the distance to the rotation axis. Poisson's equation becomes: 



Po d (6p\ d^> d£ n d^ e 
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where 



^3 {lr ( r g r (g - r 2 r (( - r 2 r gg + 2rr 2 - r 2 g r a - r 2 r g cot 6») , 



d 2 gg + coted g - 



sin 2 9 



(A. 13) 
(A. 14) 



As was pointed out in Sect. 2.2 the relative Lagrangian density perturbation, 6p/po, can be eliminated in favour of the relative 
Lagrangian pressure perturbation, 6p/P{), thanks to the adiabatic relation, Eq. ([6j. 



A.2. Mechanical boundary condition 



As explained in Sect. 2.5 the mechanical boundary condition is obtained by calculating the dot product between and Euler's 
equation, and cancelling out the vertical gradient of 6p/Po- Furthermore, the quantity is eliminated through the continuity 
equation, and the terms Q g and d g Q g vanish at the surface. In spheroidal components, one obtains: 
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Appendix B: Lagrangian perturbation to the effective gravity 



As was explained in Sect. |3.2| the Lagrangian perturbation to the effective gravity, 6g e ff, is deduced from the vectorial Lagrangian 
perturbation to the effective gravity, Sg e ff, via the relation <Sg e ff = —n ■ 6g e s, where n is the outward normal at the surface. 
Furthermore, 5g e ff includes the Lagrangian perturbation to the gradient of the gravitational potential and the acceleration of a 
particle tied to the surface, resulting from the oscillatory motions. After adding and subtracting £ ■ V \s£l 2 e s ) in order to introduce 
the equilibrium effective gravity, a vectorial expression is obtained in Eq. ( 34 1, and is reproduced here for convenience: 



Sg eS = -V¥ + g ■ Vg eff + (to + mQ.) 2 ^ - 2i(o) + mQ)Q x£-Qx(fix£)-£-V (sQ 2 e 5 ) 



(B.l) 



In what follows, we will go through the above equation one term at a time in order to obtain explicit expressions for the dot product 
of n with each one. 

The Eulerian perturbation to gravity is obtained through tensor analysis: 
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(B.2) 



where we have used the relation n 



! ^- Tn E i > . Furthermore, we have used Einstein's summation convention on repeated indices. 



Before dealing with the next term, it is useful to introduce the contravariant components of the Lagrangian displacement, which 
we distinguish from the components given in Eq. ( |A.4[ > by placing a tilde over the top: 

(B.3) 



where (E^, Eg, E<p) is given in Eq. (A.l i. We also introduce the covariant components of the effective gravity: 

geff - ~g ( t- - g g t . 



(B.4) 



From Eq. (|3T|, it is straightforward to see that g = — 3,-Po/po- Furthermore, g^f = L —Yp 

I | 1 i ( r2+r e) " 

In tensonal notation, the term n- {£ ■ Vg- e ff } becomes: 



geS and g e g = at the stellar surface. 
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?(ftgf-I*ygf)^. 



(B.5) 



where g 1 ' = E' ■ E> denotes the contravariant components of the metric tensor, and T*. = (diEj) ■ E k - - (djE k ) ■ Ej the Christoffel 
coefficients. Equation (B.5 1 can be simplified if we use the following relation: 



V • g eB = -g ij (d i8 f - T* jg f) = -47TGP0 + V ■ (stfe s ) , 



(B.6) 
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which is a modified form of Poisson's equation. The result is: 



where 



« • {£ ■ Vg eff ) = - — ~~[J2 [-^Gp + V ■ (sn 2 e s )] + ? (g!%gf - Agf) + f (-g*%gf + Bgf)} , (B.7) 
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and where we have made use of the following simplifications: gf 1 = dggf 1 — on the surface, g 1 ' = g'' and r*. = 1%. The above 
expression can then be re-expressed in terms of f r , and g e ff to yield: 
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The term V ■ \sGl 2 e s ) takes on the following expression for a general rotation profile, Q. = Q(£, 8): 



V • (s£l 2 e s ) = 2Q 2 + 



, sin 8 (r sin 8 - r a cos ff) i 



d ( {n 2 ) + sin 8 cos 8d B (p 2 ) . 



For a cylindrical rotation profile, Q = £2(s), it becomes: 

V • (s£l 2 e s ) = 2Q 2 + s0, (n 2 ) . 

The term (w + mO) 2 n • £ is given by: 



(o) + mO) 2 n • f = (<u + mo) 2 ^ ^ 
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The Coriolis term is given by: 
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(B.ll) 



(B.12) 



(B.13) 



(B.14) 



The last two terms are more conveniently treated together. They take on the following expression for a general rotation profile: 

r£ r - r g & 



-i2x(fix£)-£-V (sQ. 2 e s ) = - 



sin 0d ( (fl 2 ) + & sin 8d e (o 2 ) 



If the rotation profile is cylindrical, they become: 

-£lx(£lxf)-£-V (iQ 2 e s ) = - (sin 6% r + cos 0&) sd, (Q 2 ) e s 
Combining all of these equations together, and remembering the minus sign, yields Eq. ( [35| . 

Appendix C: Cancelling of simplified disk-integration factors 



(B.15) 



(B.16) 



Given the simplified form of the disk-integration factors given in Eq. (45 i (see Sect. |3.7| >, it turns out that some of these cancel out 
regardless of inclination. To see this, one needs to start with an explicit form for Eq. B5j: 

% j X-o f<p- 6 -f(6) Sf(8)e'" 1 ^ l+ ' a " [r (sin i sin 8 cos <p + cos i cos 8) - rg (sin i cos 8 cos <p - cos i sin 8)] r sin #d#d0 



D(i) cos(wf + i/r) 



7Tfl 2 q (ST) 



(CI) 

where i// is a suitably chosen phase, (5T) is given in Eq. (45 1, and f(8) corresponds to the visibility curve (i.e. the border between 
the visible and hidden side of the star). The function / obeys the following symmetry: f(8) + f(n - 8) — n. We have made use of 
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Eqs. ( 16 1 and (22 1 in obtaining an explicit expression for e b s . • dS. The integral in Eq. (C.l i is then split into two equal halves and 
the second half is modified according to the variable changes ff — n — 6 and <// = (f> — n: 

D(j)cos(<Uf + ^) = %{ I 6f(9)e' m "' + ' a "[r (sin i sin 9 cos + cos i cos 9) - rg (sin i cos 9 cos <p - cos i sin 9)] r sin 0d6tl0 

[Je=o J0=-/(fl) 

+ (-I)" 1+1 5f(7r-0 , )e ,m0 ' +m '[r(sin/sin0'cos0' + cos i cos 60 

J <?'=() Jip'=f(6')-27T 

- r ff (sin / cos 0' cos 0' - cos i sin 9')] r sin flWd^' | / (2nR^ <#T>) , (C.2) 

where we have made use of the relation f(9) + f(n -9) — n. The two halves can be combined to give a single integral over the entire 
stellar surface only if 5T{9) = {-\)"' +l 5T{n-9). If we assume this is the case, we can then see under conditions the integral vanishes. 

We specifically look at the integration over </j, which now is between the bounds and 2n. Remembering that L e lm ^ cos (f>d<p - 
if \m\ + 1 and that j^" e'^dcp — if m + 0, we deduce the second condition for cancelling the disk-integration factor, i.e. \m\ > 2. 



Appendix D: Normalisation of multi-colour visibilities 

In order to find a normalisation which minimises the distances between a set of multi-colour visibilities, we start with the cost 



function given in Eq. (47 1 and include an additional term so as to enforce the constraint £f =1 Wf = 1. Without loss of generality, we 



work with the normalised components, V., instead of the original ones: 
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where A represents a Lagrange multiplier. Setting the derivatives, dJ/dAj, dJ/dWi, and dJ/dA, to zero leads to the following system: 



V J ■ W 

Aj = —. r, 

V J ■ V J 

N 

(N-A)W = ^AjV 3 , 

7=1 

W W = 1, 



(D.2) 



(D.3) 



(D.4) 



where we've used vectorial notation for conciseness and where A ■ B represents the dot product 2, = i A;B;. Using the normalised 



components V. allows us to simplify Eq. (D.2 1 to Aj — V J ■ W. When combined with Eq. (D.3 1, this yields: 

N 

(N - A)W = ^ (V J ■ W) V j . 



(D.5) 



7=1 



This last equation is in fact an eigenvalue problem where N - A is the eigenvalue and W the eigenvector. In order to determine 
which eigensolution yields the minimal value for J, we develop the cost function as follows: 

N N N 

J = J] ■ V } - 2A ~jV j ■ w + w ■ w } = Z i 1 ~ A2 j\ = N ~ Yu^P (D ' 6) 

7=1 7=1 7=1 

where we've used the simplifications V J ■ V J — 1, Aj ; = V J ■ W, and WW- 1. In order to simplify the term we calculate 

the dot product between W and Eq. (|D.3|i: 



(tY - A)W ■ W = N - A = Aj (V j ■ W) = J^Aj. 



(D.7) 



7=1 



7=1 



Hence, 



J = N - (N - A) = A. (D.8) 

Therefore, the minimal value of A (and hence the maximal value of N - A) corresponds to the minimal value of J. The vector W 
is therefore the principal component of the vector set (V J ^ and can be found either via a singular value decomposition of the 
associated matrix, or more simply through a power iteration. 



